Quasi-normal mode analysis in BEC acoustic black holes 
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We perform a quasi-normal mode analysis of black hole configurations in Bose- Einstein conden- 
sates (BEC). In this analysis we use the full Bogoliubov dispersion relation, not just the hydro- 
dynamic or geometric approximation. We restrict our attention to one- dimensional flows in BEC 
with step-like discontinuities. For this case we show that in the hydrodynamic approximation 
quasi-normal modes do not exist. The full dispersion relation, however, allows the existence of 
quasi-normal modes. Remarkably, the spectrum of these modes is not discrete but continuous. 
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I. INTRODUCTION 

Analogue models of gravity [J, Q allow the study of 
certain geometrical aspects of linearized general relativ- 
ity (GR). When the analogue models are based on an 

' intrinsically quantum system, such as Bose-Einstein con- 

! densates (BECs) |, H @, this opens up the possibil- 
ity of studying some corrections to this geometric de- 

\ scription. These corrections — encoded in the form of 
non-relativistic or modified dispersion relations, similar 
to what is expected in a wide range of scenarios for grav- 
ity at Planck-scale energies, leading to possible violations 

' of Lorentz invariance [6| — are then of a quantum origin, 

, but they express themselves in a classical manner. In gen- 
eral, these corrections have to be taken into account in 

\ any treatment involving high energies (high frequencies). 
In BECs in particular, in the presence of horizons, the 

\ deviations from the geometric or hydrodynamic regime 
are important even at low energies 0]. 

In a previous article 0], the stability of several one- 
dimensional configurations in BECs was studied. These 
configurations were chosen in analogy with the geome- 
tries associated with gravitational black holes. Here we 
will elaborate further on black hole configurations, which 
were shown to be devoid of instabilities, and examine 
their quasi-normal or relaxation modes. The analysis 
we carry out is analogous to the standard quasi-normal 
mode (QNM) analysis in gravity [j| According to 
GR, when black holes are perturbed, they emit gravi- 
tational waves in the following way. The initial regime 
depends mainly on the concrete form of the perturba- 
tion. Then, an oscillatory phase is attained in which the 
precise form of the emitted wave depends only on the 
properties of the black hole itself. Finally, a polynomial 
tail characterizes the return to equilibrium. In the inter- 
mediate phase, a discrete set of complex frequencies are 
excited: the quasi-normal modes of our concern. Since 
quasi-normal modes are modes of decay, i.e. of energy 
dissipation, they are outgoing. In GR this means that 
their group velocities are directed outwards both in the 
asymptotic region and at the horizon. Moreover, as we 



will illustrate, in GR this also automatically implies that 
they are non-normalizable or divergent. 



A general investigation of quasi-normal modes in ana- 
logue black holes was carried out in [nj. As we will 
briefly point out, in standard GR in (1 + 1) dimensions, 
quasinormal modes do not exist. In higher dimensions, 
a discrete spectrum appears [1, [t| . Similar results were 
found in [1(| for analogue black holes. Here the essential 
element that we add is the influence of a non-relativistic 
contribution to the dispersion relation on these QNMs 1 . 
With the full Bogoliubov dispersion relation in BEC, 
not only does a QNM-spectrum appear even in the one- 
dimensional case studied here, but this spectrum turns 
out to be continuous. This result suggests the straight- 
forward speculation that in (3 + 1) dimensions, for both 
analogue and GR black holes, there will also be continu- 
ous regions of QNMs when modifications of the dispersion 
relation with respect to the hydrodynamic or relativistic 
case are taken into account. 



The structure of this paper is as follows. In the next 
two sections we will review the steps that lead from the 
general description of a BEC to the simplified formal- 
ism for one-dimensional configurations with abrupt dis- 
continuities, and highlight the modified dispersion rela- 
tion. In Sec. IIVI we will discuss the boundary conditions 
for our problem, and compare with the standard quasi- 
normal mode analysis of gravitational black holes. Then, 
in Sec. [V] we present and discuss our main results. 



Modifications of the quasi-normal mode spectrum of black holes 
due to Lorent violations have also been investigated in 
but from a different perspective. They introduce spontaneous 
Lorentz-symmetry breaking terms in an effective Lagrangian via 
non-zero expectation values of different tensor fields in the vac- 
uum, thereby restricting attention to a modified Dirac equation 
in a Schwarzschild background. However, these terms do not 
lead to modified dispersion relations of the type analyzed here. 
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II. PRELIMINARIES 

In this section we will briefly review the description of 
background BEC configurations and their perturbations, 
with particular emphasis on the modified dispersion re- 
lation. The main aim here is to fix notation. General re- 
views of BECs can be found e.g. in [l2|,[l3[ while the type 
of configurations discussed here were amply described in 

A dilute gas of interacting bosons can be described in 
terms of quantum field operators \I/(x) [and ^(x)] that 
annihilate [and create] particles. When this gas conden- 
sates, \& can be separated into a macroscopic wave func- 
tion ip and a field operator Lp describing quantum fluc- 
tuations: 4" = ip + (p. The macroscopic wave function tp 
satisfies the Gross-Pitaevskii (GP) equation 

m iW x) = v 2 + Kxt(x) + g \ip\A i>{t, x), 



dt 



2m 



(2.1) 



where m is the boson mass, V cx t the external potential 
and g a coupling constant which is related to the corre- 
sponding scattering length a through g = Ant^a/m. 

This equation can be expressed in terms of hydrody- 
namic quantities such as the local speed of sound c and 
the velocity of the fluid flow v. This proceeds by first 
introducing the Madelung representation for the order 
parameter in the GP equation: 



4, 



(2.2) 



Here n is the condensate density, fi the chemical potential 
and 9 a phase factor. Next, as we are interested in the 
acoustic perturbations of the BEC, we linearize in n and 
9: 

n(x, t) = n (x) + g _1 ni(x, t), (2.3a) 
0(x,t) = o (x)+0 1 (x,t). (2.3b) 

Then, the following set of equations for the linear per- 
turbations is obtained from the complex GP equation 



d t ni = -V • (niv + c 2 V6>i) 



d t 9i 



1 9 

-v • V6»i -m + -fV 



c 2 V 



(2.4a) 
(2.4b) 



which arc equivalent to the Bogoliubov equations, and 
where £ = — is the healing length, and c and v are 



c 2 — gno/m, v — V9o/m. 



(2.5) 



The last term on the right-hand side in (I2.4b[) is the so- 
called "quantum potential" term. 

At this point, let us briefly recall that in the hydrody- 
namic limit, which is obtained by neglecting the quantum 
potential term just mentioned, an effective Lorentzian 
metric can be constructed of the type: 



-v 
1 



(2.6) 



which historically formed the onset of the interest in hy- 
drodynamic analogue models for gravity (for an overview, 
see [l| and references therein). 

From now on, consider one-dimensional profiles. The 
equations (|2.4p can be solved in a homogeneous region (c 
and v constant) by writing the perturbations in plane- 
wave form 

ni(x,t) = Ae i( - kx - ut \ (2.7a) 
9 x {x,t) = Be l{kx ~^\ (2.7b) 

This immediately leads to the local dispersion relation 2 



(w - vkf = c 2 k 2 + \c 2 ek\ 



(2.8) 



This dispersion relation is quartic in k. It should be 
compared to the hydrodynamic dispersion relation that 
is obtained by neglecting the quantum potential: 



{lu ~ vk) 2 = c 2 k 2 



(2.9) 



This hydrodynamic dispersion relation is of the usual 
quadratic form. 

However, here we are not so much interested in this 
hydrodynamic limit, but precisely in deviations from it. 
Therefore we will always work with the full dispersion re- 
lation, and only use the hydrodynamic limit as a means 
of comparison. The key point is that, because the full 
dispersion relation is quartic, when writing the mode cor- 
responding to a certain frequency lu: 



,-twt 

3 



A, 



(2.10) 



there are now four contributions to this u^-mode, stem- 
ming from the four values of k associated to each value 
of lu. In the hydrodynamic limit there are obviously only 
two contributions to each w^-mode: 



hi, 



: Lu/(v ± c). 



(2.11) 



In the configurations analyzed in this paper, the addi- 
tional, "non-hydrodynamic" modes will turn out to play 
a crucial role even at low frequencies. 



III. BACKGROUND CONFIGURATIONS AND 
NUMERICAL METHOD 

We are looking for quasi-normal modes of black hole 
configurations in one-dimensional BECs of the type rep- 
resented in fig. [1] These profiles are piecewise homoge- 
neous with a single discontinuity. In order to find these 



2 The term 'dispersion relation' is sometimes reserved strictly for 
an equation of the form u> = io(k), while an equation such as ||2.8|I 
would then have to be called e.g. the 'characteristic equation'. 
For the sake of clarity, we will always use the term 'dispersion 
relation' to refer to an equation of the form 1 12. 811 or 112.91 . 
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FIG. 1: (Color online) Flow and sound velocity profile consist- 
ing of two homogeneous regions with a step-like discontinuity 
at x = simulating a black hole-like configuration in a BEC. 
The solid blue line represents the speed of sound c, the dashed 
red line the fluid velocity v. In the upper part of the picture, 
the negative value of v indicates that the fluid is left-moving. 
For x > 0, the fluid is subsonic since c > \v\. At x < it has 
become supersonic. At x = 0, there is a sonic horizon. 



relaxation modes, first, we will need a way to describe 
the mode solutions that characterize the profiles under 
study. The idealized form of these profiles allows a sim- 
plified description in each homogeneous region. Then, 
the various sections have to be linked through match- 
ing conditions at the discontinuity. After finding all the 
mode solutions, the quasi-normal modes will be those 
that moreover satisfy the adequate boundary conditions. 
We will describe the matching conditions in the current 
section, together with a straightforward algorithm for nu- 
merical implementation. The boundary conditions will 
be the topic of the next section. 

Integrating (|2.4[) in an infinitesimal neighbourhood of 
the discontinuity (say x = 0) and simplifying with the 
use of the continuity equation {vc 2 = constant) leads, for 
our choice of a piecewise homogeneous background, to 
the following set of matching conditions Q : 



\e x ] = o, [ c 2 aA] = o, 



ni 

r-2 



0. 



[d x ni 



0, 



(3.1a) 
(3.1b) 



where the brackets indicate, for example, 

[0\] = Q\\ x =a+ - 0i| x =o-- 

In each homogeneous section, we look for a solution for 
a particular frequency u> as a sum of £;-modes of the form 



3 



vk 



ic 2 k 2 



j_ e i(kjx—uit) 



(3.2) 
(3-3) 



where the sum is over the four fcy-values associated to u> 
in that homogeneous region. 



The conditions (|3.ip can be written in matrix form 

(3.4) 



AyAj = 0. 



Since there are eight free parameters Aj for each dis- 
continuity (four in each region), and only four condi- 
tions, matching is in principle always possible. However, 
there are additional conditions: the boundary conditions 
which we will discuss in the next section. These boundary 
conditions can be added to the matrix Ay, to obtain a 
(4 + N) x 8 matrix Ajj , with N the number of constraints 
on the coefficients Aj that follow from the boundary con- 
ditions. Depending on the value of N, we can have the 
following cases. 

• If N < 4, then there will always exist a non-trivial 
solution Aj, since there are more degrees of free- 
dom than constraints. Therefore, a region of the 
complex w-plane where N < 4 will represent a con- 
tinuous region of eigenfrequencies. 

• If N = 4, then we have a 8 x 8 system of equa- 
tions and there is a non-trivial solution Aj only if 

det(A) = 0. Generally speaking, this means that 
in a region where N — 4, we might at the very 
most expect isolated discrete eigenfrequencies, but 
no continuous zones of solutions. 

• If N > 4 then the system can be split into two (or 
more) 8x8 subsystems, each with a subdetermi- 
nant Aj. For a non- trivial solution to exist, each of 
these subsystems should fulfill the previous condi- 
tion det(Aj) = 0. 

To summarize, we can define a non-negative function 
F(u>) such that 

• F(u) = if N < 4, 

• F(u) = |det(A)| if AT = 4, 

• F(u) = E |det(Ai)| if N > 4. 

The eigenfrequencies of the system are those values of u> 
for which F(u>) = 0, and since the quasi-normal modes 
are relaxation modes they will be located in the lower 
half complex u> plane. 



IV. QUASI-NORMAL MODE ANALYSIS AND 
BOUNDARY CONDITIONS 

The essential ingredient that is still missing to calcu- 
late the quasi-normal modes in these BEC black hole 
configurations are the boundary conditions. We will first 
take a look at the boundary conditions that are used to 
calculate QNMs in GR. This can be immediately extrap- 
olated to BECs in the hydrodynamic limit. However, as 
we will show, in this hydrodynamic limit in (1 + 1) di- 
mensions, there exist no QNMs. Finally, we will come to 
the boundary conditions with the full dispersion relation. 
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A. Quasinormal modes of a gravitational black hole 

Quasi-normal modes are simply the modes of energy 
dissipation of a perturbed object. Gravitational black 
hole QNMs in particular are the relaxation modes that 
characterize the pulsations of the black hole after pertur- 
bations initiated internally in the spacetime surrounding 
the black hole. 

The fact that these modes decay in time means that 
they can be represented by complex frequencies u> with 
Im(w) < 0. The fact that they are relaxation modes 
means that they will have to be outgoing in the asymp- 
totic region. At the horizon, QNMs are also required to 
be outgoing because in GR nothing can escape through 
the horizon (note that in QNM terminology "outgoing" 
means directed towards the exterior of the spacetime re- 
gion connected to an asymptotic observer, in other words 
at the horizon this means directed towards the singular- 
ity) . So the boundary conditions that are imposed in GR 
to find the QNMs of a black hole are simply that these 
modes should be outgoing both in the asymptotic region 
and at the horizon. 

Since these QNMs are outgoing, they are also diver- 
gent. Indeed, because of the quadratic dispersion rela- 
tion, say J 1 = c 2 k 2 , we have for the group velocity v g 
and for the imaginary part of k: 

v g — ±c, Im(fc) = lm(uj)/v g 

and therefore, since QNMs have lm(w) < 0: 

sign[Im(fc)] = -sign[u fl ]. 

To sum up, QNMs in a gravitational black hole are found 
by imposing the boundary condition that they should be 
outgoing, and because of the quadratic dispersion rela- 
tion they are also automatically divergent. 



B. BEC hydrodynamic limit: absence of 
quasi-normal modes in one-dimensional flows 

In BECs in the hydrodynamic limit, decaying outgoing 
modes are also automatically divergent, again because of 
the quadratic dispersion relation. Indeed, 

V g = v ± c, Im(fc) = Im(uj)/vg 

and therefore sign[Im(fe)] = — sign[u g ] as before. Let us 
see whether QNMs can exist in this case. 

In the hydrodynamic limit, the quantum potential 
term in (|2.4[) is neglected, and so upon integration of the 
remaining terms in these equations, only two matching 
conditions remain: 

[0i] = O, [vn 1 +c 2 d x 9 1 }=0. (4.1) 

Note that this second equation cannot be simplified as in 
(|3.1[) , because this requires the use of a third matching 



condition which came precisely from integration of the 
quantum potential. 

The hydrodynamic solution for a particular 
frequency ui, 

ni = Aie^'-^ + A 2 e l{k2X ~ u>t \ (4.2) 
0i = B l e i{klX ^ t) + B 2 e i{k2X - uit \ (4.3) 

has to satisfy Aj = i(u> — vkj)Bj. Assume that 
the configuration is completely subsonic. Then, for 
x < 0, ki = u//(v + c) = fc ; ~ is ingoing, whereas 
&2 = u)/(v — c) = k~ ut is outgoing, and vice versa for 
x > 0, where we write k^ and k^ ut respectively. In terms 
of the associated coefficients Bj , the matching conditions 
then become 

B+ ut + B+=Br+B- at , (4.4) 
c+ (B+ Ut - B+) = c- (Br n - B- Ut ), (4.5) 

or in matrix form 

For a solution to be outgoing in both regions, both B^ n 
and Bj~ should vanish. This means that the coefficient 
c + + c~ connecting B^ with B~ ut must be zero. But 
since c is always a positive integer, this is impossible. 
Therefore, there are no quasi-normal solutions in the hy- 
drodynamic limit. 

The above calculation is strictly valid only for com- 
pletely subsonic profiles. When the left-hand side is su- 
personic, the hydrodynamic approximation breaks down 
in the passage through the sonic horizon [7J • Both modes 
in the supersonic region become outgoing at x — > — oo, 
while there is still one outgoing mode in the right asymp- 
totic (subsonic) region. Since the matrix in (|4. 6[) is well 
defined, it seems then that some QNMs exist for these 
configurations. However, this is a spurious effect caused 
by the subsonic to supersonic step-like discontinuity in 
combination with the hydrodynamic regime. Precisely 
because of the presence of the horizon, we know that 
none of the modes at the left-hand side can be connected 
to the outgoing mode at the right-hand side. So QNMs 
cannot exist. From a point of view more similar to GR, 
when the (x < 0) region is supersonic, it becomes discon- 
nected from the subsonic region (x > 0). But then, just 
like in GR, the requirement that QNMs must be outgoing 
in the left asymptotic region should be replaced by the 
requirement of being outgoing at the horizon. Because 
of the homogeneity of the profile, none of the modes can 
be outgoing both at the right asymptotic region and at 
the horizon. So in any case, QNMs cannot exist in the 
hydrodynamic limit in (1 + 1) dimensions, regardless of 
the sub- or supersonic character of the profile. 

Note that in GR in (1 + 1) dimensions, a more gen- 
eral argument prohibits the existence of QNMs. Indeed, 
all two-dimensional metrics are conformally flat. In ad- 
dition, the dAlembertian equation in two dimensions is 
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conformally invariant [14(. Therefore, all the solutions 
of the system are conformally equivalent to plane waves, 
and so they don't satisfy the requirement of being out- 
going both in the asymptotic region and at the horizon. 
Hence QNMs cannot exist in (1 + 1) GR. The situation is 
completely different in (3 + 1) dimensions, in which GR 
black holes exhibit a discrete QNM spectrum. 
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C. Boundary conditions with the full dispersion 
relation 

As we just discussed, in BECs in the hydrodynamic 
limit, decaying outgoing modes are automatically diver- 
gent. But when looking at the full dispersion law, this 
relation is no longer valid. Note that for our choice of 
profiles, we have two asymptotic regions, and moreover 
no prohibition of crossing the horizon in any direction, 
since the permeability of the horizon is an essential fea- 
ture of modified dispersion relations [7]. An outgoing 
mode is then one that has a group velocity v g < for 
x < and v g > for x > 0, where v g is obtained from 
the full dispersion relation: 

Clearly, there is no immediate connection with the sign 
of Im(fc) anymore, so the automatic divergent character 
of outgoing modes is lost. 

In any case, quasi-normal modes should be outgoing, 
and this remains the essential boundary condition. Ad- 
ditionally, one might wonder whether the QNMs will be 
convergent, or divergent after all (in spite of the modifi- 
cation of the dispersion relation). An easy way to check 
this is by applying convergence as an additional bound- 
ary condition and comparing the resulting spectrum with 
the one obtained without this additional condition. 

To sum up, quasi-normal modes should be outgoing. 
In GR and in the hydrodynamic approximation, this au- 
tomatically implies a divergent behaviour in the asymp- 
totic regions. However, in both these cases, there are 
no QNMs in (1 + 1) dimensions (as opposed to (3 + 1) 
GR black holes with their discrete QNM spectrum). In 
BECs with the full modified dispersion relation, outgo- 
ing modes are no longer automatically divergent. The 
boundary condition to determine QNMs in the (1 + 1) 
BEC black hole configurations under study is therefore 
that they should be outgoing. Nonetheless, we will check 
whether they are convergent or divergent. 



V. RESULTS AND DISCUSSION 

In this section we present and discuss results for 
the quasi-normal mode spectrum of black holes in one- 
dimensional BECs obtained with the full dispersion rela- 
tions. 



FIG. 2: Plots illustrating the quasi-normal mode spectrum 
for a one- dimensional black hole configuration in BEC. The 
left-hand part of the picture represents the number N of con- 
straints following from the boundary condition that QNMs 
should be outgoing. The right-hand part shows the function 
log[.F(w)] defined in Sec. IIHI where F(u) = corresponds 
to quasi-normal modes. The QNM spectrum consists of the 
continuous region where N < 4. [The numerical values used 
for these plots, in units such that the healing length £ = 1, 
are c = 1; V = 0.7 in the subsonic region and v = 1.8 in the 
supersonic region.] 



Fig|5] illustrates the QNM spectrum for an acoustic 
black hole containing a single discontinuity, subsonic at 
the right-hand side and supersonic at the left. Remember 
from pj that such a black hole did not possess unstable 
eigenfrequencies, so it indeed makes sense to look for sta- 
ble modes. 

As discussed in the previous section, the boundary con- 
dition that determines the existence of QNMs is that they 
should be outgoing in both asymptotic regions. Surpris- 
ingly, as can be seen from Fig[21 not only does such a 
QNM spectrum indeed show up, but this spectrum is 
actually continuous. 

The next question is whether these modes are con- 
vergent or not. In Fig [31 we have imposed convergence, 
respectively at both sides, in the right asymptotic region 
only, and in the left asymptotic region (the singularity) 
only. The QNM spectrum disappears as soon as the con- 
vergence condition is imposed at either side. In other 
words, the quasi-normal modes are divergent. Remem- 
ber from the previous section that for outgoing hydro- 
dynamic modes, this divergent character is automatic. 
However, we have also seen that when only retaining the 
hydrodynamic modes, there are no QNMs. So although 
the QNMs that we have just found are divergent, the in- 
clusion of the non-hydrodynamic modes is essential for 
their existence, and so their divergent character is not 
"automatic" , i.e. it is not due to a strict relation between 
the signs of v g and Im(fc). As a matter of fact, we have 
checked that the QNMs disappear when the convergent 
fc-contributions are excluded. 

We have also calculated the QNM spectrum for other 
black-hole like configurations that were shown to be de- 
void of instabilities in and we have checked that a 
similar continuous spectrum of QNMs appeared in all of 
them. This indicates that the continuous character of the 
QNM spectrum is an essential consequence of the mod- 



G 



N(<fl) -log [F(»)] 




Re((jl) He(ffl) 



FIG. 3: Alteration of QNM spectrum resulting from imposing 
the additional boundary conditions of (from top to bottom) 
convergence in both asymptotic regions, in the right asymp- 
totic region only, and in the left asymptotic region only. Since 
F(oj) never reaches zero, it is as if the QNM spectrum of Fig. [2] 
had disappeared. In other words, the QNMs are divergent in 
both asymptotic regions. [Same numerical values for c and v 
as in Fig. [1] 

ified dispersion relation, and not just of the particular 
characteristics of a precise configuration. 

To sum up, there is a continuous region of 
frequencies ui such that they have associated modes 

4 
3 = 1 

which are outgoing at both ends. These modes are di- 
vergent, just like in the hydrodynamic case, although the 
fact that here non-hydrodynamic modes come into play 
means that this relation cannot be derived directly from 
the dispersion relation. 

VI. CONCLUSIONS 

We have examined the quasi-normal modes of one- 
dimensional black hole configurations in BECs, and laid 



particular emphasis on the importance of the deviation 
of the full dispersion relation with respect to a relativis- 
tic or acoustic one. The full dispersion relation is quartic 
in BECs, as compared to the usual second order of its 
hydrodynamic and general relativistic equivalent. With 
such quadratic dispersion relations, outgoing relaxation 
modes are automatically divergent. In particular, this 
argument is valid for QNMs in GR black holes. But the 
modification of the dispersion relation means that this 
relation is no longer automatic. Nevertheless, the quasi- 
normal modes in the systems discussed here are also di- 
vergent. 

More importantly, the QNM spectrum that was ob- 
tained in these systems with a modified dispersion rela- 
tion consists of a continuous region in the complex fre- 
quency plane. The importance of this result is easy to 
see when taking into account that both in GR and in 
BECs in the hydrodynamic limit, QNMs simply do not 
exist in (1 + 1) dimensions, while in higher dimensions the 
quasi-normal modes of a GR black hole form a discrete 
spectrum. So due to the modification of the dispersion re- 
lation, the QNM spectrum in one-dimensional BEC black 
holes changes from non-existent to a continuous region of 
frequencies. 

We would like to end this paper with the straightfor- 
ward speculation that the discrete spectrum which is ob- 
tained in the usual quasi-normal mode analysis in stan- 
dard (3 + 1) GR, will also develop continuous regions 
when taking the modification of the dispersion relation 
at high energies into account. The existence of these con- 
tinuous bands could become a signal of trans-Planckian 
physics encoded in the emitted spectrum of gravitational 
waves. 
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